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To test a theory of cosmic microwave background fluctuations, it is natural to expand an anisotropy 
map in an uncorrelated basis of linear combinations of pixel amplitudes — statistically-independent 
for both the noise and the signal. These S/A-eigenmodes are indispensible for rapid Bayesian 
analyses of anisotropy experiments, applied here to the recently-released two-year COBE dmr maps 
and the firs map. A 2-parameter model with an overall band-power and a spectral tilt i<at describes 
well inflation-based theories. The band-powers for all the dmr 53, 90, 31 a-\-b GHz and firs 170 GHz 
maps agree, {(1.1 ± 0.1) x 10 - } ' , and are largely independent of tilt and degree of (sharp) S/N- 
filtering. Further, after optimal S/A-filtering, the dmr maps reveal the same tilt-independent large 
scale features and correlation function. The unfiltered dmr 53 a-\-b index i<at + 1 is 1.4 ± 0.4; 
increasing the S/A-filtering gives a broad region at (1.0-1.2)±0.5, a jump to (1.4-1.6)±0.5, then a 
drop to 0.8, the higher values clearly seen to be driven by S/A-power spectrum data points that do 
not fit single-tilt models. These indices are nicely compatible with inflation values (~0.8-1.2), but 
not overwhelmingly so. 
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PACS NOs: 98.80.Cq, 98.80.Es, 98.70.Vc 

The importance of the large-angle COBE dmr [1,2] 
and the balloon-borne firs [3] detections for testing the- 
ories of cosmic structure formation can hardly be over- 
stated. For theories in which the multipole components 
of the radiation pattern at m are Gaussian-distributed and 
statistically-independent, as in inflation-based models or 
via the central limit theorem, we only need to deter- 
mine the power spectrum, Ct = £(£ + 1) (|a£ m | 2 )/(27r). 
A phenomenology characterized by just two parameters, 
a broad-band power (fit) dmr and a broad-band tilt ^at 
is an excellent approximation for a large class of models: 
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defines the "logarithmic integral" of a function ft- For 
the dmr and firs maps, the broad-band filter W ' t encodes 
beam-smearing and effects from pixelization [9]; for other 
experiments it also encodes switching information. The 
virtue of the band-power [4-6] (fit)yj- is that it is insen- 



sitive to the specific shape of Ct 
a Qrms ps" varies considerably: 



whereas the oft-used [1] 
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a=0.55 for firs. The broad-band tilt v&t measures de- 
viation from scale-invariance in AT: Ut={£ + 2)" AT (1 + 
0((£ + i) -2 )- It is usual to use n s = v&t + 1 to charac- 
terize the slope [1], because for pure scalar metric fluctu- 
ations in a cold dark matter Universe with almost no 
baryons, n s is the primordial (e.g., post-inflation) in- 
dex [7]. However, the standard 5% baryon content gives 
i>^ T = 0.15 for an initially untilted primordial spectrum 
[8,5]; and primordial tilts, usually negative, can arise 



from gravitational wave as well as density fluctuation 
modes, but with n s > 0.7 required to get the observed 
cluster abundance (e.g., [8,6]). 

The observations (AT /T) p from the p th pixel of an 
Npix CMB anisotropy experiment is given to us in terms 
of an average value (AT/T) and a correlation ma- 
trix C'opp' , with possibly off-diagonal components as 
well as the diagonal <t|j . The theory we are test- 
ing is characterized by a correlation matrix Ct pp > (here 
=X[CtWtPt(cos6ppi)], where 6 pp i is the angle between 
pixels p and p'). If signals are not Gaussian-distributed, 
higher order correlation tensors are required, and such 
theories can only be well-analyzed by simulation. 

For this exploration of the dmr maps, I used: one lower 
resolution scale than the original maps, 5.2° c.f. 2.6°; a 
Galactic latitude cut |6| > 25° (leaving 928 pixels); dipole 
and average subtractions were done after cuts, but be- 
fore resolution-lowering; the (non-Gaussian) revision of 
the dmr beam [9], with corrections for digitization and 
pixelization; a correction linear in the off-diagonal Co 
components to test sensitivity to residual correlation be- 
cause COBE actually measures a AT between pixels 60° 
apart. For the 168 GHz firs map, 2.6° rather than 1.3° 
pixels were used (but I checked that the two give the 
same answer), leaving 1070 pixels; a 3.9° beam was used, 
which includes pixelization corrections. 

A full Bayesian analysis of maps requires frequent in- 
version and determinant evaluations of N p i x x N p i x cor- 
relation matrices, the sum of all C'xpp' in the theoreti- 
cal modelling plus the pixel-pixel observational error ma- 
trix Copp'- To facilitate this, I expand the pixel values 
(AT/T)p into a basis of "signal-to-noise" eigenmodes for 
the maps in which the transformed noise and transformed 
(wanted) theoretical signal we are testing for do not have 
mode-mode correlations, i.e., are orthogonal. This can 
always be done, no matter what the experiment [4]. Com- 
plications are associated with the removal of averages, 



dipoles, etc. and the existence of secondary signals in 
the data, both of which do couple the modes. A model 
for the various contributions that make up the observed 
data is then 
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where R is a rotation matrix. The modes are sorted 
in order of decreasing Sy TV-eigenvalues, £TR,k, so low k- 
modes probe the theory in question best. With uniform 
weighting and all-sky coverage, the SyTV-modes are just 
the independent Re(ai m ) and Im(ai m ), with the lowest 
I having the highest £TR,k, hence k ~ (£ + l) 2 . With 
inhomogeneous pixel coverage, Galactic cuts, dipole sub- 
tractions, etc., they are more complex. 

This expansion is a complete (unfiltered) representa- 
tion of the map. The sum of £| over bands in S 1 / TV-space 
defines a S/N-power spectrum which gives a valuable 
picture of the data; an example is shown in Fig.l. The 
maps can have an arbitrary average, dipole and possibly 
quadrupole (because the Galactic contribution may con- 
taminate the signal's quadrupole). This is modelled by 
Ck, with 4 or 9 components, which I take to be Gaussian- 
distributed with very wide width, i.e., with a uniform 
prior probability. In S/N-space, the variance JCkk' has 
off-diagonal terms affecting the largest scale S 1 / TV-modes 
and complicating the interpretation of Fig.l for small k, 
but in the Bayesian analysis these unknowns can imme- 
diately be marginalized (by integrating over all cj;). 

Fig.l also shows anomalies relative to the power-law 
theory at high k (see [5] for more examples). I model 
these residuals by an excess pixel noise with an amplifi- 
cation factor r. For firs, r ~ 0.25. For the 2-year data, 
the most likely r < 1% for the recommended dmr noise 
level for all maps, often zero for z^at > 0.2. For dmr, I 
include it anyway because of the high sensitivity of the 
full Bayesian method, but the results differ very little 
if I set r = 0. Features in Fig.l that do not fit the 
smooth theoretical model (around k ~ 200 for 53a+b) 
represent unknown components, denoted by res^, which 
could point to a better model power spectrum to try, or 
residuals left over from the various dmr map-processing 
techniques. In general, the res^ covariance 1Z would be 
nondiagonal. The nature of the low and high S 1 / TV-modes 
has been probed with angular correlation functions [6]. 
These reveal the expected result that high 5'/TV-modes 
do indeed involve collective large scale pixel combina- 
tions, while low SyTV-modes involve destructive interfer- 
ences from nearby pixels that are insensitive to large scale 
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FIG. 1. The (5/TV) band-power (mean of £ k — 1 over a band 
in k) for the dmr 2-year 53a+b maps, using the (a) i/at=0 and 
(b) j/at=0.5 modes. The observational error bars (smaller 
than the points for k > 30) include a large scale power cor- 
rection. The <0.01 values at high k are actually slightly neg- 
ative. The small solid points joined by the central dotted 
curve give band-averages of £rR,k, with the upper and lower 
dotted curves giving theoretical variances for this binning; all 
three scale with (CA , , thus can be raised or lowered to best 
fit the data. The horizontal line of solid triangles shows the 
most likely residual noise (tiny for (a), off-scale for (b)). The 
Weiner-filter is the long-dashed curve. The 90a+b (open tri- 
angle) and 31a+b (x) data points shown are uniformly scaled 
until their £tr,u agrees with 53a+b (which they do): thus, 
given the larger error bars, the agreement in power among 
maps is reasonably good. 



structure in the maps, but are quite sensitive to physics 
inside the beam, whether from systematic effects or true 
'white' noise on the sky. Thus S / N -modes form an ideal 
set for filtering. 

The first step in the Bayesian method is the construc- 

1/2 

tion of a joint likelihood function in (Ci) dmr , v&t and 
r (adopting a uniform prior probability in each of these 
variables). Integrating over r (marginalizing it) allows 

-I/O 

one to construct v^t~{Ci) dmr contour maps. Marginal- 
ization over ^at or cuts along a given ^at (e.</., 0) gives 
a {Ct)J mr distribution; marginalization over {Ct)J mr gives 
a i/^T distribution. The Table gives 50% Bayesian prob- 

1/9 

ability values for (Ci) dmr and vat , with 'one-sigma' error 
bars from the 84% and 16% Bayesian values, 'two-sigma' 
from the 97.5% and 2.5% values [11]. The {Ci) l J^ r for all 
of the 53, 90 and 31 GHz dmr maps and firs map agree 
at better than the one-sigma level, along fixed i/^T lines, 
among different z^at, and when i/^T is marginalized, 
essentially independently of the degree of S/ TV-filtering 
(cuts at 400 and 200 are shown). There is also no dis- 
cernable power in the dmra—b maps. These band-powers 
compare with: the original first year dmr number, us- 



ing correlation functions (CF) [1], 0.97±.28; a later up- 
date [12], 0.97±.16 using the Boughn-Cottingham single 
quadratic statistic (for ^at=0); the two-year data, using 
CF [2], 1.00±.10 (for v AT =0) and 1.02±;|f (for the most 
probable v^t)', and the firs CF result [15], 1.08±0.3. 

The tilts shown in the Table (with one and two sigma 
errors) are not determined with the same precision, nor 
can one attach the same confidence to the values, be- 
cause of the variation with filtering. In particular, the 
high point for the 53a+b S 1 / TV-power in Fig.l just be- 
low k = 200 (which comes from 53a and not from 53b) 
drives the index up somewhat, as do the two slightly high 
points around 100 (which 90a+b has also), while the low 
point at 70 drives it down. Thus the most probable in- 
dex is sensitive to filtering, as Fig.l makes clear; it also 
shows a high index power law is not a particularly good 
description of the excess power. Another significant fea- 
ture is that the low quadrupole drives the index higher 
(the '+q' results). I do not actually have to subtract 
the quadrupole since the large uniform prior allows the 
analysis to attribute what it cannot fit to c k . I agree 
with the dmr team [14] that the lower indices obtained 
by assuming Galactic interference in the quadrupole are 
the better numbers to adopt. Decreasing the beam-size 
also lowers ^at a bit [6]. The ^at results compare with: 
the original first year dmr value [1], 1.2_ 07 ; the two-year 
CF result [2] 1.6_ 06 ; power spectrum (PS) values using 
quadratic estimators [13] for one-year, l-7_ ' 5 , and two- 
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(for 3 < £ < 19), 1.31^ (for 3 < £ < 30), 
and using Bayesian methods with linear multipole fil- 
tering [14], l.ltoi, (for 3 < £ < 30). For the first 
year 53a+b map, I get 1.8_ y^ i with no quadrupole, 
2.0_oY 1 *o wl th, because the k ~ 200 structure is more 
pronounced [6]. The firs CF value [15] is 1-0_ 5 includ- 
ing the quadrupole, whereas I get l.S^o'g'.j/g for firs,+q, 
and rest; is large enough that direct 5'/#-filtering can 
give even higher values [6]. 

I have found sharp 5'/#-filtering preferable for sta- 
tistical analysis [16], but smooth preferable for clean- 
ing noise to look for robust map-independent features. 
Fig. 2 shows optimally-filtered maps and their tempera- 
ture autocorrelation functions, C{9). The maps are equal 
area projections [17] about the pole. Floating averages, 
dipoles and quadrupoles can be added to bring them into 
better agreement. This Weiner-filtering is an immediate 
byproduct of the 5'/#-eigenmode expansion [5,18]: given 
observations £& , the mean value and variance matrix of 
the desired signal s k are [19] 
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The mean field (sj;|£) is the maximum entropy solution. 
The operator multiplying £& is the Weiner filter. The 
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FIG. 2. The first column shows unfiltered 100° diameter 
dmr maps centered on the North Galactic Pole, the second 
column shows them after Weiner-filtering, the third the South 
Pole version, assuming a i/at=0 C^-spectrum, with the nth 
contour as noted, positive contours heavier than negative ones 
and the zero contour heaviest of all. The x-axis points to the 
Galactic center, the j/-axis to qcjNP x qac, the Supergalactic 
plane is long-dashed, and the ecliptic plane is short-dashed; 
all are uncorrelated with the features. The second row en- 
tries are: the first-year Weiner-filtered GNP result; two-year 
GNP, GSP results, but with i<at = 1. The maps have been 
smoothed by a 3° Gaussian filter. In the lower figure, correla- 
tion functions for 53a+b GHz (open circles) and 90a+b GHz 
(solid triangles) are contrasted with the Weiner filter correla- 
tion functions for both cases, each for i<at = 0, 1. All but the 
unfiltered 90a+b GHz one agree. 
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fluctuation, Asj, = s fc — (sfc|<f), of the signal about the 
mean is realized by multiplying a vector of N p i x indepen- 
dent Gaussian random numbers by the square root of the 
variance matrix. 

The Weiner filter depends upon the (Ci) dmr and ^at 
we choose: I use the maximum likelihood value for 
(Ci) dmr ; and both maps and correlation functions are rel- 
atively insensitive to ^at for large scale features. When 
the noise is large, as it is for these maps, it is the higher 
£ power that is preferentially removed by optimal filter- 
ing. Thus the theoretical fluctuation Asj, would have to 
be added to the Weiner-filtered maps of Fig. 2 to give a 
realistic picture of the sky given the data and the theory. 
That is the maps are too smooth, 31 GHz the most, 53 
GHz the least. The optimally-filtered correlation func- 
tions for both 53 and 90 GHz converge to the raw 53 GHz 
correlation function, and are ^AT-independent (but this 
does not address how statistically significant such a C{9) 
is for the given z^at)- Not shown are the C{9) for the op- 
timally filtered dmr a— b difference maps, but these are 
nicely zero. Also not shown, but very gratifying, is that 
the main features are evident in both A and B channel 
maps, including in the noisy 90a and 31a, b. 

Although correlation functions are quadratic combina- 
tions of pixel amplitudes and thus more complicated than 
the linear (Gaussian-distributed) statistics used here and 
in [14], they are effective filters of small-angle high- 
frequency systematic effects: a pixel-error enhancement 
(r) only contributes to the zero angle bin of C{9), and 
the low S / N -modes at high k have a range below about 
9 jwhm , the beamsize [6]. Thus ^at determined with C{9) 
should be a reasonable estimate if bins within 9f w hm are 
not included (although this increases the error bars). 

This letter shows that the overall anisotropy power 
COBE has discovered exists in all channels and is very 
robust. What drives the index ^at high is more indica- 
tive of leftover residuals than suggestive of high-tilt pri- 
mary anisotropies which are problematic for inflation- 
based models. Although more computer-intensive sta- 
tistical exploration is still needed to determine the best 
vat, especially when the four full years of dmr data are 
released, it seems unlikely the errors will ever be small 
enough that a strong argument against inflation based 
solely on dmr can be made. This research was supported 
by NSERC and the Canadian Institute for Advanced Re- 
search, and owes much to the nice job the dmr team did 
with the two-year data set. 
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